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ABSTRACT 

In  a  review  of  studies  on  the  seismic  phase  Lg,  we  describe  its  particle 
motion,  dispersion,  spectral  content,  mode  of  propagation,  and  magnitude-scale; 
we  also  tabulate  the  regional  velocity,  attenuation,  and  propagation  efficiency 


for  this  seismic  phase. 


The  characteristics  of  Lg-wave  propagation  in  the  eastern  United  Spates 
are  compared  with  those  in  different  regions  of  the  Soviet  Union.  Possible 
discriminants  such  as  (i)  Lg  vs.  P  amplitudes,  (ii)  Lg/P  amplitude  ratios  as 
a  function  of  distance,  and  (iii)  Lg  energy  ratios  are  found,  similar  to  atten¬ 
uation  and  group  velocity,  to  be  highly  dependent  on  the  propagation  path.  The 
valid  application  of  these  quantities  to  the  problem  of  earthquake-explosion  dis¬ 
crimination  will  therefore  require  regional  studies  more  detailed  than  previously 


assumed. 


A  re-evaluation  of  the  magnitude-yield  relation  and  an  examination  of 
physical  parameters  which  may  be  relevant  to  the  estimated  yield  of  underground 
nuclear  explosions  were  performed.  The  preliminary  results  indicate  that  (i) 
the  m  vs.  yield  relation  shows  regional  differences  and  dependence  on  the  source 
medium,  and  (ii)  the  collapse  volume  and  the  diameter  of  the  collapsed  crater  are 
usually  proportional  to  the  estimated  yield. 
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INTRODUCTION 

The  content  of  this  report  is  divided  into  three  topics:  (i)  a  review 
of  the  available  studies  on  the  seismic  phase  Lg,  (ii)  a  comparison  of  regional 
wave  propagation  in  the  US  and  the  USSR,  and  (iii)  a  preliminary  re-evaluation 
of  the  magnitude-yield  relation  and  an  examination  of  the  physical  parameters 
which  may  be  relevant  to  the  estimated  yield  of  underground  nuclear  explosions. 

Beginning  with  the  current  contract  period,  we  will  be  reporting  a  series 
of  reviews  on  (i)  the  seismic  phases  (Lg,  Rg,  Pg,  Pn,  Sn)  that  are  potentially 
useful  to  the  discrimination  of  explosions  from  earthquakes  at  regional  distances, 
and  (ii)  the  spectral  characteristics  of  underground  nuclear  explosions.  A  lack 
of  critical  reviews  on  these  crucial  subjects  has  motivated  us  to  undertake  this 
ambitious  project.  We  hope  to  achieve  three  goals  through  the  reviews:  (i)  to 
compile  and  categorize  the  available  observations  into  accessible  format,  (ii) 
to  summarize  the  theoretical  development  in  an  overview  fashion,  and  (iii)  to 
emphasize  the  features  that  are  related  to  the  problems  of  earthquake-explosion 
discrimination.  In  this  report,  we  will  present  a  review  on  the  seismic  phase 
Lg.  The  review  is  subdivided  into  7  topics:  (A)  particle  motion  and  dispersion, 
(B)  regional  velocity,  (C)  spectral  content,  (D)  wave  guide  and  mode  of  propaga¬ 
tion,  (E)  attenuation  and  propagation  efficiency,  (F)  magnitude-scale  based  on 
Lg,  and  (G)  others  (Sn-to-Lg  conversion,  application  to  the  earthquake-explosion 
discrimination  problem,  and  search  for  oceanic  Lg). 

A  comparative  study  of  regional  wave  propagation  in  the  eastern  United 
States  and  different  regions  of  the  Soviet  Union  is  presented  in  the  second  part 
of  this  Semi-Annual  Technical  Report.  Five  topics  were  selected  to  assess  the 
feasibility  of  directly  comparing  the  characteristics  of  regional  seismic  waves 
in  the  US  and  the  USSR,  and  to  evaluate  their  relative  importance  to  the  problem 
of  earthquake-explosion  discrimination.  The  topics  are:  (i)  Lg  vs.  P  amplitudes, 
(ii)  Lg /P  amplitude  ratios  as  a  function  of  distance,  (iii)  Lg  group  velocity, 

(iv)  Lg  energy  ratios,  and  (v)  Lg  attenuation. 

In  studying  regional  seismic  wave  propagation,  we  often  encounter  the 
problem  of  how  to  calibrate  a  magnitude-yield  relation  at  regional  distances. 

This  problem,  although  quite  fundamental  in  nature,  is  by  no  means  an  easy  one 
because  a  well-determined  magnitude-yield  relation  requires  a  clear  knowledge  of 
(i)  the  source  size,  (ii)  the  amplitudes  of  seismic  waves  at  different  distances, 
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(iii)  the  effects  of  crustal  structure  at  the  source  and  the  receiver,  and 

(iv)  the  effects  of  the  propagation  path.  The  last  part  of  this  report  re¬ 
examines  this  relation;  it  also  describes  the  preliminary  results  from  analyses 
of  several  physical  parameters  that  are  related  to  the  yield  of  underground 
nuclear  explosions. 


The  purpose  of  this  review  is  threefold:  (i)  to  provide  a  summary  of 
the  available  observations  on  Lg,  (ii)  to  present  the  theoretical  developments 
in  an  overview  fashion,  and  (iii)  to  clarify  or  comment  on  what  appears  to  us 
to  be  confusing  concerning  the  interpretation  of  Lg. 

The  name  Lg  was  assigned  by  Press  and  Ewing  (1952)  in  their  pioneering 

study  on  this  seismic  phase.  "L"  because  the  particle  motion  was  predominantly 

of  Love  or  transverse  type,  and  "g"  because  the  wave  was  believed  to  propagate 

in  the  granitic  layer  of  the  crust,  and  was  therefore  considered  a  surface-wave 

counterpart  of  the  near-earthquake  body  waves  Pg  and  Sg.  These  authors  sum¬ 
marized  the  properties  of  Lg  (for  propagation  paths  in  North  America)  succinctly 
in  the  abstract  of  their  1952  paper: 

"Surface  shear  waves  (Lg)  with  initial  period  1/2  to  6  seconds  with  sharp 
commencements  and  amplitudes  larger  than  any  conventional  phase  have  been  recorded 
for  continental  paths  at  distances  up  to  6000  km.  These  waves  have  a  group  velo¬ 
city  of  3.51  ±  0.07  km/sec.  and  for  distances  greater  than  20°  they  have  reverse 
dispersion.  For  distances  less  than  about  10°  the  periods  shorten  and  Lg  merges 
into  the  recognized  near-earthquake  phase  Sg." 

This  and  later  investigations  of  Lg  also  point  out  that  (i)  the  wave  is  not  ob¬ 
served  after  approximately  100  km  of  propagation  in  the  oceanic  crust,  (ii)  the 
particle  motion  may  contain  a  substantial  amount  of  longitudinal  and  vertical 
components,  and  (iii)  the  observations  may  be  explained  by  a  collection  of  Airy 
phases  of  hiqher  mode  Love  and  Rayleigh  waves. 

The  terms  of  Sg  and  Lg  were  used  to  refer  to  different  waves  in  some  earlier 
studies.  Although  both  terms  referred  to  high-frequency  shear  waves  in  the  con¬ 
tinental  crust,  the  distinctions  were  based  on  differences  in  the  observed  fre¬ 
quency  content,  the  distances  of  observation,  and  the  interpretation  in  their 
mode  of  propagation.  Sg,  which  is  analogous  to  its  compressional-wave  counterpart 
Pg,  referred  to  the  direct  shear  arrival  at  short  epicenter  distances;  while 
Lg  referred  to  the  superposition  of  normal  modes,  with  frequencies  slightly  lower 
than  those  of  Sg,  at  epicentral  distances  greater  than  about  10°  (Press  and 
Ewing,  1952).  {^There  has  been  considerable  confusion  concerning  the  definitions 
of  Pg  and  Sg.  These  terms  replaced  the  P  and  S  of  Mohorovicic  (191*1)  for 
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typographical  convenience  (pg.  86  of  Jeffreys,  1976)  and  the  supposed  association 
with  the  granitic  layer  of  the  crust.  While  the  definition  of  P  referred  to  the 
direct  compressional  arrival  at  short  distances  with  a  velocity  of  about  5.5  km/ 
sec  (cf.  Fig.  18-1  of  Richter,  1958),  the  original  data  was  obtained  at  distances 
over  150  km.  Explosion  data  from  California  indicated  that  direct  compressional 
arrivals  at  120  km  within  the  epicenter  had  a  velocity  near  6.34  km/sec.  ~he 
Californian  researchers  consequently  suggested  the  notation  "p"  for  the  direct 
wave  at  short  distances  and  "P"  for  the  compressional  wave  with  a  velocity 
around  5.5  km/sec  (p.  286-287  of  Richter,  1958).  The  consensus  at  the  present 
seems  to  be  the  use  of  the  nomenclature  P  for  direct  compressional  waves  and  the 
terms  "Pn"  and  "Pg"  for  occasions  when  two  distinct  arrivals  with  velocities 
around  8. 0-8. 4  km/sec  and  5. 4-5. 7  km/sec  are  observed.')]  In  view  of  the  con¬ 
sensus  on  the  terminology  of  P-,  Pg-,  and  Pn-  waves  and  the  arbitrary  distinction 
between  Sg  and  Lg,  we  are  in  favor  of  calling  the  direct  shear  arrival  "S"  and 
reserving  the  term  "Lg"  for  shear  waves  with  group  velocities  around  3.5  km/sec 
at  epicentral  distances  where  Sn  (or  the  mantle-refracted  S)  becomes  the  first 
shear  arrival.  In  this  report,  the  term  "Lg"  will  refer  to  both  the  "Lg"  and 
the  "Sg"  cited  in  earlier  seismological  literature.  In  the  following  sections, 
we  will  attempt  to  summarize  and  discuss  previous  studies  on  the  observations  and 
interpretations  of  the  Lg  phase.  We  have  divided  the  literature  available  to  us 
into  7  topics:  (A)  particle  motion  and  dispersion,  (B)  regional  velocity,  (C) 
spectral  content,  (D)  wave  guide  and  mode  of  excitation,  (E)  attenuation  and 
propagation  efficiency,  (F)  magnitude-scale  based  on  Lg,  and  (G)  others. 

(A)  Particle  motion  and  dispersion. 

Press  and  Ewing  (1952)  describe  the  particle  motion  of  Lg  in  the  following 

words : 

"...During  the  first  cycles  the  waves  have  approximately  equal  amplitudes 
on  all  three  components,  but  the  transverse  horizontal  rapidly  gains  amplitude 
and  becomes  several  times  larger  than  the  other  two  within  about  30  seconds.  Ap¬ 
proximately  1  minute  after  the  commencement  of  the  phase,  the  amplitude  on  the 
transverse  component,  having  reached  a  value  many  times  larger  than  that  of  S 
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or  SS  on  any  component,  begins  to  decrease  gradually,  but  does  not  drop  to  a 
value  comparable  with  that  of  SS  until  about  30  minutes  later,  the  period  then 
being  of  the  order  10-14  seconds.  The  group  velocity  for  the  latter  part  of 
this  phase  is  certainly  less  than  about  2  km/sec,  the  lower  limit  being  uncertain 
As  for  Eurasian  events  recorded  at  Uppsala  and  Kiruna,  Bath  (1954) 
reports  that  the  particle  motion  of  Lg  was  primarily  transverse  and  was  often 
observed  at  two  different  group  velocity  windows:  Lg-j  ,  at  3.54  ±  0.06  km/sec  and 
L 92  at  3.37  +  0.04  km/sec.  Lehmann  (1953)  states  that  there  was  "considerable" 
vertical  motion  involved.  All  the  authors  mentioned  above  agreed  that  both 
the  horizontal  and  the  vertical  components  of  particle  motion  were  present  in  the 
Lg  phase.  Herrin  and  Richmond  (1960)  used  a  ray-approach  analysis  to  explain 
the  particle  motion  of  Lg.  Their  calculations  indicate  that  a  strong  SV  type 
motion  (i.e.  with  longitudinal  and  vertical  components  of  motion)  would  be  present 
with  the  SH-type  motic  initially;  but  during  the  later  part  of  the  wave  train 
where  the  angle  of  incidence  for  the  rays  presumably  becomes  less  steep,  energy 
leakage  to  the  bottom  layers  due  to  SV-to-P  conversion  would  occur  and  the  SV- 
motion  tends  to  decrease  faster  than  that  of  the  SH-motion.  The  results  of  this 
analysis  are  in  agreement  with  the  observations  of  Oliver  et  al .  (1955),  but 
do  not  agree  with  their  own  observations  at  Dallas  for  earthquakes  in  southwestern 
United  States  and  Mexico  where  strong  SV-motion  continued  throughout  the  Lg  wave- 
train.  Herrin  and  Richmond  also  estimated  the  partitioning  of  energy  between 
SV  and  P  waves  at  different  angles  of  incidence;  Herrin  (1961)  pointed  out  some 
errors  in  their  partitioning  of  energy  and  corrected  them.  By  correlating  the 
vertical  component  to  the  longitudinal  component  of  the  Lg  particle  motion, 

Sutton  et  al.  (1967)  found  out  that  the  particle  motion  of  Lg  from  underground 
nuclear  explosions  and  small  earthquakes  tended  to  be  either  transverse  or  mixed. 

Aside  from  the  qualitative  comparison  of  Press  and  Ewing  between  the  vertical 
and  horizontal  components  of  displacement,  there  are  several  other  reports  on  their 
»  relative  amplitudes.  For  the  Lg  amplitudes  generated  by  the  nuclear  explosion  GNOME 

in  a  salt  mine  of  New  Mexico,  Romney  et  al.  (1962)  note  that  the  displacements 
cn  all  three  components  were  approximately  equal.  But  for  earthquakes  in  the 
northeastern  U.S.  -  southeastern  Canada  regions  recorded  at  North  American  stations, 
t  Street  (1976)  reports  that  the  maximum  sustained  horizontal  component  of  Lg 
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consistently  exceeded  the  vertical  component  by  a  factor  of  3.  For  all  epicentral 
distances  in  Iran,  the  resultant  horizontal  motion  of  Lg  at  1  sec  was  usually 
twice  that  of  the  vertical  component  (Nuttli,  1980a).  Bath  (1956),  however, 
found  some  Lg  waves  with  no  vertical  particle  motion  at  all. 

Although  Press  and  Ewing  (1952)  suggested  the  possibility  of  using  higher 
mode  surface  waves  to  interpret  the  Lg  phase,  Oliver  and  Ewing  (1957)  were  the 
first  to  calcualte  the  dispersion  curves  of  higher  mode  Rayleigh  waves  and  use  them 
to  explain  the  longitudinal  and  vertical  components  of  Lg  particle  motion.  In 
a  later  paper,  Oliver  and  Ewing  (1958)  computed  the  dispersion  curves  from  simple 
earth  models  for  higher  mode  Love  and  Rayleigh  waves  and  found  that  the  M^-mode 
(1st  shear  mode)  and  the  second  Love  mode  had  similar  velocities  at  the  same  period, 
which  may  explain  the  simultaneous  arrivals  of  the  vertical,  longitudinal,  and 
transverse  components  of  ground  motion  for  Lg.  Dispersion  curves  and  particle 
motions  of  higher  mode  Love  and  Rayleigh  waves  were  computed  for  realistic  earth 
models  by  Brune  and  Dorman  (1963),  and  later  including  the  effects  of  sphericity 
into  the  earth  models  by  Kovach  and  Anderson  (1964).  Brune  and  Dorman  also  com¬ 
puted  synthetic  seismograms  for  the  transverse  component  of  Lg.  The  results  of 
these  authors  confirm  the  hypothesis  of  Oliver  and  Ewing.  Knopoff  et  al .  (1973) 
presented  further  evidence  to  identify  the  transverse  component  of  Lg  motion  as 
higher  mode  Love  waves  by  (i)  computing  the  relative  spectral  excitations  for 
double-couple  sources  at  different  depths,  and  (ii),  constructing  synthetic 
seismograms  for  the  higher  mode  Rayleigh  waves  and  identified  them  as  the  long¬ 
itudinal  and  vertical  components  of  Lg  motion. 

The  particle  motion  of  the  1st  shear  mode  (M?)  was  computed  by  Oliver  and  Ewing 
(1957)  to  be  retrograde  elliptical;  the  same  authors  later  reported  that  obser¬ 
vations  from  an  Arctic avent  (5/25/1950,  8:34:32;  65.5°N,  151. 5°W)  recorded  at 
Palisades,  confirmed  their  previous  theoretical  results  on  the  particle  motion 
(Oliver  and  Ewing,  1958).  Barley  (1978)  traced  the  particle  motion  of  nigner 
mode  Rayleigh  waves  (2.0sec  T  3.5sec)  for  the  group  velocity  wndow  3.0 

to  3.5  km/sec,  and  found  it  to  be  retrograde  elliptical.  This  result  was  predicted 
by  the  theoretical  calculations  of  Panz3  et  al .  (1972)  for  the  first  three  higner 
Rayleigh  modes;  these  authors  also  found  that  at  a  given  period  the  ellioticity 
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(defined  as  the  ratio  of  the  longitudinal  component  of  particle  motion  to  the 
vertical  component)  increased  with  decreasing  mode  number.  For  a  shield 
structure  with  a  low  velocity  channel  (LVC)  in  the  upper  mantle,  they  found 
that  at  periods  less  than  4  sec.  the  ellipticity  for  the  third  higher  Rayleigh 
mode  was  greater  or  equal  to  0.7,  whereas  the  ellipticity  for  the  fundamental 
and  the  first  two  higher  Rayleigh  modes  was  greater  or  equal  to  1.0. 


(B)  Regional  velocity. 

Table  I  is  a  summary  of  Lg  velocities  which  were  published  in  journals 
and  reports  available  to  us.  Whenever  possible,  we  tried  to  include  infor¬ 
mation  pertaining  to  the  measurements  of  the  velocity,  such  as  the  location 
of  the  seismic  events  and  recording  stations,  the  type  of  instrument  used  to 
record  the  events  (horizontal  or  vertical  component,  short  or  long  period,  etc.), 
and  the  period  of  the  Lg  waves  at  which  the  measurement  was  made.  Although 
the  majority  of  the  references  cited  did  not  specify  their  method  of  measure¬ 
ment.  we  deduced  from  their  figures  that  most  reported  velocities  were  measured 
at  the  initial  stage  of  the  coda  when  a  visible  change  in  wave  frequency  or 
amplitude  could  be  observed,  either  on  the  long-  or  short-period  instruments. 

The  measurements  of  Pomeroy  and  Nowak  (1978),  however,  were  made  at  the  amplitude 
maxima  of  the  Lg  coda  which  seemed  to  be  more  unstable.  Differences  in  the 
method  of  measurement  and  the  recording  instrument  may  account  for  the  apparent 
discrepancy  between  the  various  reports.  While  measurements  at  the  beginning 
of  the  coda  probably  correspond  to  the  Airy  phase(s)  of  higher  mode  surface  waves 
with  the  fastest  group  velocity,  measurements  at  the  amplitude  maxima  probably 
coincide  with  the  group  velocity  window  where  several  Airy  phases  overlap. 

Whereas  the  former  is  indicative  of  the  average  properties  of  the  wave  guide, 
the  latter  which  tends  to  be  slower  than  the  former,  is  probably  not  only  more 
diagnostic  of  the  detailed  structure  of  the  wave  guide  but  also  informative 
concerning  the  relative  excitation  of  the  various  modes  at  the  source  (Knopoff 
et  al .  ,  1974).  We  would  like  to  explore  this  possible  aspect  of  Lg  in  a  future 
study. 
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(C)  Spectral  Content. 

The  only  sources  known  to  us  on  the  spectral  content  of  Lg  are  derived 
from  Street  et  al.  (1976)  and  the  Soviet  seismological  literatures  (e.g.  Anto¬ 
nova  et  al . ,  1978;  Nurmagambetov ,  1974).  The  studies  on  Lg  propagation  in 
the  USSR  were  compiled  and  summarized  in  a  report  by  Shishkevish  (1979). 

Street  et  al .  derived  their  data  from  over  300  short-period,  vertical 

component  recordings  of  78  earthquakes  in  the  central  U.S.  In  the  period 

range  they  analyzed  (approximately  0.05  -  10  sec.),  the  amplitude  spectra  gen- 

_2 

erally  indicate  a  falloff  of  oj  between  the  flat  portions  at  the  long-  and 
short-  period  ends.  Their  spectra  were  corrected  for  the  effects  of  instru¬ 
ment  response,  but  not  for  the  anelastic  attenuation  of  the  path. 

The  frequncy  selection  seismograph  stations  (ChISS)  of  the  USSR  have  enabled 
the  spectral  analysis  of  Lg  to  become  a  routine  procedure.  Their  results,  com¬ 
monly  plotted  as  log  (A/T)  vs.  log  (1/T),  generally  display  peaks  at  short 
epicentral  distances.  The  peak  is  shifted  towards  lower  frequencies  as  epicentral 
distance  increases.  This  dependency  of  spectral  peak  on  epicentral  distance 
is  also  a  function  of  propagation  path.  In  these  studies,  the  frequency  ranged 
from  0.3  to  approximately  20  Hz  while  the  epicentral  distance  spanned  from  30 
to  3000  km.  The  falloff  in  their  velocity  amplitude  spectra  (i.e.  displacement 
amplitude  spectra  multiplied  by  frequency)  is  also  dependent  on  epicentral 
distances:  at  epicentral  distances  around  350  km,  the  falloff  ranges  from 
slightly  greater  than  one  to  approximately  two;  whereas  at  epicentral  distances 
greater  than  about  1000  km,  the  falloff  remains  less  than  3.  Since  these 
measurements  of  Lg  spectral  content  did  not  take  the  effects  of  geometrical 
spreading  and  anelastic  attenuation  into  account,  the  spectral  characteri sti cs 
measured  at  short  epicentral  distances  were  probably  more  representative  of  the 
source  spectra  and  a  spectral  falloff  of  about  2  could  be  taken  as  representative 
of  the  source  falloff  for  the  displacement  amplitude  spectra  of  Lg  waves.  The 
high-frequency  spectral  peaks  observed  in  the  USSR  is  probably  an  artifact  of 
the  velocity  spectra  plot;  that  is,  the  spectral  peak  will  disappear  if  the  plot 
is  converted  into  a  displacement  amplitude  spectra. 
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(D)  Wave  guide  and  mode  of  excitation. 

Press  and  Ewing  (1952)  are,  again,  the  first  ones  to  point  out  that 
"...  Lg  is  a  wave  which  is  confined  to  a  surface  or  near-surface  layer  by  wave- 
guid^action. . . "  based  on  the  observed  velocity  and  large  amplitudes.  Subsequent 
theoretical  studies  tend  to  support  their  claim  although  this  conclusion  is  not 
reached  without  its  share  of  confusion.  In  a  study  of  Lg  waves  in  Eurasia, 

B4th  (1954)  observed  a  correlation  between  hypocentral  depth  and  the  energies 
contained  in  Lg1  and  Lg2-  That  is,  the  energy  of  Lg1  generally  decreased 
with  increasing  hypocentral  depth,  whereas  the  energy  for  Lg2  reached  a  maximum 
when  the  source  depth  was  around  45  km.  He  attributes  the  difference  in  energy 
distribution  to  several  crustal  channels  or  layers  which  transmitted  waves  at 
different  group  velocities.  This  claim,  although  sound  when  interpreted  in 
terms  of  Airy  phases  with  different  group  velocities,  led  to  two  unexpected 
results  when  viewed  from  the  perspective  of  channel  waves.  Firstly,  terminologies 
for  waves  which  supposedly  propagated  in  different  channels  of  the  crust  and 
upper  mantle  proliferated  (e.g.  Bath,  1958).  Secondly,  several  low-velocity 
channels  in  the  crust  and  upper  mantle  came  to  be  used  as  explanations  for  the 
efficient  propagation  of  the  various  channel  waves  (Gutenberg,  1955;  Bath,  1956, 
1958). 

Based  on  the  dispersion  curves  of  higher  mode  Love  and  Rayleigh  waves, 

Oliver  and  Ewing  (1957,  1958),  Brune  and  Dorman  (1963),  and  Kovach  and  Anderson 
(1964)  found  it  possible  to  explain  the  frequency  content  and  the  group  velocity 
of  Lg  waves  by  using  the  Airy  phases  of  the  higher  modes.  Kovach  and  Anderson 
(1964)  also  point  out  that  the  modes  observed  "...depend  on  the  period  range 
being  studied  and  the  depth  of  the  source..."  and  that  variations  in  the 
velocity  and  period  of  the  observed  Lg  depended  on  the  positions  of  the  Airy 
phase,  wnich  in  turn  depended  on  the  elastic  parameters  of  the  propagation 
path.  If  the  interpretation  of  Lg  waves  as  superpositions  of  higher  mode 
surface  waves  is  correct,  then  we  would  expect  an  additional  dependence  on 
the  source  radiation  pattern.  At  periods  greater  or  equal  to  5  sec.,  radiation 
patterns  of  the  first  higher  Love  and  Rayleigh  modes  compare  favorably  with 
calculated  results  (Mitchell,  1973  a,b).  The  observations  of  Sutton  et  al . 

(1967)  on  short-period  (0. 5-2.0  sec)  Lg  waves,  however,  indicate  that  "... 


there  seems  to  be  no  systematic  difference  in  the  short-period  energy 
radiation  pattern  between  the  underground  nuclear  explosions  and  the 
earthquakes  ..."  and  that  the  pattern  of  the  energy-contours  (or  con¬ 
tours  based  on  the  maximum  amplitude)  could  be  better  explained  by  a 
correlation  with  the  major  tectonic  provinces  of  the  United  States. 

Since  the  modal  composition  of  Lg  at  short  periods  is  a  combination  of 
many  higher  modes,  the  observed  amplitudes  may  not  be  diagnostic  of  the 
radiation  pattern  of  the  individual  modes.  Also,  scattering  is  prob¬ 
ably  more  important  for  short-period  waves  and  its  effects  more  likely 
to  mask  any  azimuthal  pattern  that  may  be  present. 

Panza  et  al.  (1972)  showed  that  the  collection  of  higher  mode  Ray¬ 
leigh  waves  could  be  separated  into  a  family  of  crustal  waves  and  a  family 
of  channel  waves  in  a  structure  containing  even  a  slight  low-velocity 
channel  (LVC)  in  the  upper  mantle.  As  it  is  implied  by  the  name,  channel 
waves  have  most  of  the  energy  in  the  LVC  and  have  essentially  zero 
energy  at  the  surface.  Crustal  waves,  on  the  other  hand,  have  most  of 
their  energy  in  the  crust;  consequently,  only  the  fundamental  mode  and  the 
crustal  waves  need  to  be  considered  for  the  excitation  of  Rayleigh  waves. 
Knopoff  et  al.  (1973)  demonstrated  that  higher  mode  Love  waves  could 
similarly  be  divided  into  crustal  waves  and  channel  waves.  For  a  structure 
without  any  LVC,  the  whole  suite  of  higher  mode  Love  and  Rayleigh  waves 
has  to  be  taken  into  account  for  the  ground  motion  of  the  Lg  waves. 

Knopoff  et  al.  (1974)  further  establish  that  the  group  velocity 

and  the  periods  of  the  Lg  stationary  phase  could  be  diagnostic  for  the 

crustal  thickness  and  the  shear  velocity  in  the  crust  and  the  upper  mantle. 

In  general,  as  the  crustal  thickness  increased,  both  the  group  velocity 

of  the  late-arriving  La  stationary  phases,  U  .  ,  and  the  period  at  U  ■  , 

3  min  min 

T  .  ,  tended  to  increase.  Increasing  tne  crustal  velocity  while  keeping 
all  other  parameters  constant  would  tend  to  decrease  T  ^n,  but  increase 
Umin,  the  magnitude  of  Lg-excitation,  and  the  general  period-content  of 
the  Lg  waves.  These  authors  also  demonstrate  that  (i)  for  thicknesses 
of  the  upper  mantle  lid  greater  than  20-25  km,  Lg  is  insensitive  to  changes 
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in  its  thickness,  and  (ii)  Lg  is  insensitive  to  the  velocity  in  the  upper 
mantle  LVC.  Panza  and  Calcagnile  (1975)  point  out  that  higher  mode  con¬ 
tribution  becomes  more  significant  as  the  period  decreases  and/or  as  the 
hypocentral  depth  decreases. 

As  for  the  low-velocity  channel  in  the  crust  and/or  upper  mantle, 

Oliver  and  Ewing  (1958)  concluded  that  it  was  not  necessary  to  explain  the 
characteristics  of  the  Lg  phase.  Knopoff  et  al.  (1973)  and  Panza  and 
Calcagnile  (1975),  based  on  more  modes  extending  to  shorter  periods,  reached 
the  same  conclusion  concerning  the  Love-  and  Rayleigh-type  motions  of  the 
Lg  phase,  respectively. 

Most  of  the  investigators  mentioned  in  this  section  would  probably 
maintain  that  the  characteristics  of  Lg  can  be  explained  by  the  anelastic 
attenuation  of  the  crust-mantle  layers,  the  frequency  response  of  the  seis¬ 
mograph  system,  and  the  superposition  of  higher  mode  surface  waves.  Ruzaikin 
et  al.  (1977),  on  the  other  hand,  state  that  they  "...remain  unconvinced 
that  normal  modes  will  allow  useful  interpretation  of  Lg  when  more  detailed 
data  on  its  structure  are  obtained  ..."  and  suggest  that  lateral  heterogeneity 
had  a  key  role  in  shaping  the  characteristics  of  the  observed  Lg.  Their 
argument  was  based  on  the  discrepancy  between  calculations  from  higher 
mode  surface  waves  which  predicted  the  duration  of  Lg  to  be  confined  in 
the  group  velocity  windows  of  approximately  3. 5-3.1  km/sec,  and  observations 
of  the  Lg  phase  which  indicated  that  its  amplitude  was  significant  in  the 
group  velocity  window  3.5  -  2.8  km/sec.  Oceanic  Rayleigh  waves  of  the 
fundamental  mode  (T  i  12  sec)  also  exhibit  similar  "stretching"  in  duration. 
These  waves  have  nevertheless  been  instrumental  in  shaping  our  present 
understanding  concerning  the  oceanic  structure.  Thus,  while  we  share  the 
belief  with  Ruzaikin  et  al.  that  heterogeneties  in  the  propagation  path  are 
important  i n  shaping  the  waveform  of  Lg,  we  also  believe  that  the  normal 
mode  theory,  when  supplemented  with  theories  or  methods  which  can  take  hetero- 
genety  in  the  path  into  consideration  (e.g.  the  scattering  theory  of  Aki  , 
1969),  will  serve  to  improve  the  explanation  for  the  Lg  phase. 
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This  section  deals  with  the  measurement  of  amplitude-diminution  as 
a  function  of  epicentral  distance;  the  title  of  the  section  reflects,  re¬ 
spectively,  the  quantitative  and  qualitative  aspects  of  it.  The  former 
refers  to  the  rate  of  anelastic  absorption  of  the  wave's  kinetic  energy 
per  unit  distance,  while  the  latter  provides  a  descriptive  measure  for  the 
efficiency  of  the  medium  in  transmitting  Lg  waves. 

In  seismological  literature,  attenuation  is  usually  measured  in  terms 
of  the  attenuation  coefficient,  T  ,  or  the  attenuation  quality  factor,  0  . 
These  two  quantities  can  be  related  via  the  following  equation: 

■n.  f 


where  f  and  u  are  the  frequency  and  the  velocity  of  the  wave,  respectively. 
For  Lg  waves,  measurements  of  and  Q  ,  compiled  in  Table  II,  have  been 
obtained  by  three  approaches:  (i)  time-domain,  (ii)  frequency-domain,  and 
( i  i i )  coda . 

The  time  domain  approach  entails  three  steps:  (i)  measure  the  wave 
amplitude  at  different  epicentral  distances,  (ii)  correct  the  amplitudes 
for  the  effect  of  geometrical  spreading,  and  (iii)  estimate  the  T  or  Q 
that  would  explain  the  falloff  of  the  amplitude  in  relation  to  distance. 
Nutt! i  (1975,  1978,  1980  a,  b)  and  Street  (1976)  chose  to  combine  steps 
(iii)  and  (ii)  together,  and  compared  the  observed  amplitudes  directly 
with  curves  that  include  the  effects  of  geometrical  spreading  and  dif¬ 
ferent  degrees  of  attenuation.  The  frequency-domain  approach  has  the 
advantage  of  being  able  to  take  the  source  radiation  pattern  into  account. 
The  procedure  used  by  Mitchell  and  coworkers,  who  have  been  the  primary 
advocates  of  this  approach  on  higher  mode  surface  waves,  is  similar  to  that 
employed  for  the  study  of  the  fundamental  mode  (Tsai  and  Aki ,  1969). 

Again,  three  steps  are  involved  in  this  procedure:  (i)  determine  the 
amplitude  apectra  for  the  fundamental  and  higher  mode  surface  waves  by 
applying  a  frequency-veloci ty  filter  (e.g.  the  multiple-filter  technique 
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of  Dziewonski  et  al.,  1969),  (ii)  estimate  a  fault-plane  solution  from  body- 
and/or  surface-wave  data,  and  (iii)  calculate  the  attenuation  coefficient  that 
would  produce  the  best  fit  between  the  observed  amplitudes  and  the  radiation 
pattern  computed  at  each  period.  To  date,  this  approach  has  been  limited 
to  the  analysis  of  the  fundamental  and  the  1st  higher  mode  (Mitchell,  1973 
a,b;  Cheng  and  Mitchell,  1980).  The  coda  approach,  which  was  derived  from 
the  scattering  theory  of  surface  waves  (Aki,  1969),  has  been  applied  suc¬ 
cessfully  to  data  from  narrow-band  seismographs  to  establish  (i)  scaling 
laws  for  local  earthquakes,  and  (ii)  estimates  of  regional  Q  (Aki  and  Chouet, 
1975;  Chouet  et.al.,  1978;  Rautian  and  Khalturin,  1978).  Herrmann  and  co¬ 
workers  recently  modified  this  method  for  data  derived  from  broadband  seis¬ 
mographs.  They  estimated  the  regional  Q  from  Lg  waves  by  measuring  (i)  the 
predominant  frequency  in  the  coda  as  a  function  of  time,  and  (ii)  the  coda 
shape  (Herrmann ,  1980,  Singh  and  Herrmann,  1979). 

The  propagation  efficiency  of  a  region  is  usually  estimated  by  measuring 
the  frequency  content  and  wave  amplitude  (usually  in  relation  to  the  level 
of  the  ambient  noise  or  the  amplitude  of  another  phase);  in  general,  three 
terms;  clear,  weak,  and  none,  are  used  to  describe  the  amplitude  of  the 
Lg  phase.  "Clear"  usually  refers  to  an  impulsive,  large-amplitude,  high- 
frequency  arrival;  "weak"  refers  to  a  drawn-out,  small,  low-frequency  arrival; 
and  "none"  is  indicative  of  completely  inefficient  Lg  propagation.  Although 
different  authors  have  set  their  standards  for  clear  and  weak  Lg  somewhat 
differently,  their  conclusions  concerning  the  propagation  efficiency  of  a 
given  region  are,  surprisingly,  quite  uniform.  A  list  of  regional  studies 
on  the  propagation  efficiency  of  Lg  is  compiled  in  Table  III. 

In  interpreting  the  inefficient  propagation  of  Lg  in  the  Tibetan 
plateau,  Ruzaikin,  et  al.  (1977)  proposed  two  explanations  which  are  probably 
applicable  to  most  areas  with  major  tectonic  boundaries.  Firstly,  a  dis¬ 
ruption,  termination,  or  vertical  displacement  of  wave  guide  (which  is  either 
the  entire  crust  or  part  of  it)  will  seriously  affect  the  propagation 
efficiency  of  Lg  waves;  secondly,  high  attenuation  in  the  crust  will  also 
be  able  to  affect  the  ability  to  transmit  Lg.  The  ocean-continent  boundary 
is  probably  a  disruption  or  termination  of  the  wave  guide  for  Lg;  disap- 


pearance  of  Lg  waves  after  crossing  approximately  100  km  of  oceanic  struc¬ 
ture  is  a  well  documented  observation  (e.g.  Press  and  Ewing,  1952;  B&th,  1954; 
etc.).  This  peculiar  property  of  Lg  waves  to  propagate  only  in  the  continental 
crust  was  used  by  Oliver  et  al.  (1955)  to  map  the  continental  structure  in 
the  Arctic  regions. 

Bcith  (1956)  and  Gutenberg  (1955)  report  that  the  Lg  phase  was  weakened 
or  disappeared  when  crossing  recent  mountain  chains.  Shishkevish  (1979),  in 
his  compilation  of  studies  on  Lg  propagation  in  the  Soviet  Union,  a^so  notes 
that  the  Lg  phase  was  attenuated  when  crossing  Tien  Shan,  Pamir-Hindu  Kush, 
and  the  Himalayas.  He  also  points  out  that  "...the  propagation  of  Lg  across 
the  Tien  Shan  is  less  efficient  when  paths  are  more  oblique  to  the  trend 
of  the  range  than  when  they  are  perpendicular  to  it  ...".  Uniformity  of  the 
structure  (Chinn  et  al.,  1980)  and  the  complexity  of  geology  (Street,  1976) 
in  the  propagation  path  are  also  considered  important  in  determining  the 
attenuation  of  the  Lg  amplitude.  In  summary,  the  presence  of  a  uniform, 
high-Q  wave  guide  is  essential  for  the  efficient  propagation  of  Lg;  in  the 
case  of  a  non-uniform  or  low-Q  wave  guide,  the  degree  of  non-uniformity  of 
the  wave  guide  and  the  length  of  propagation  in  it  are  both  important  in 
determining  the  fraction  of  Lg-energy  that  will  be  observed. 


(F)  Magnitude-scale  based  on  Lg. 

Since  Lg  is  often  found  to  be  the  largest  phase  at  regional  distances, 
it  ' s  natural  that  a  magnitude-scale  based  on  Lg  amplitude  would  become 
important  to  studies  on  regional  seimsicity.  Based  on  LRSM  reports  from  78 
underground  nuclear  explosions.  Baker  (1970)  proposed  a  general  formula  of 
the  form. 


%  =  log10  (A/T)  +  Q(T’  A  }  +  S(T)  (2) 
to  calculate  the  magnitude-scale  from  Lg  amplitudes.  Q  (T,a)  represents  a 
correction  term  for  the  attenuation,  and  S(T)  is  a  term  for  station  correction. 
Baker  obtained  an  expression  for  Q(T,a),  as  a  sixth  degree  polynomial  of 


distance,  by  minimizing  the  difference  between  log^A/T)  and  the  reported 
m^  for  each  event;  he  also  assigned  tentative  corrections  for  each  station, 
calculated  by  Baker  indicates  less  scatter  than  the  reported  . 

Nuttli  (1973)  formulated  a  magnitude  scale  for  Lg  while  studying  its 
attenuation  in  the  eastern  United  States.  He  assumed  that  the  term  Q(T,a) 
in  equation  (2)  has  the  form  C(T,A)  log^A,  and  subsequently  found  two 
magnitude  formulae,  applicable  at  different  distance  ranges,  for  1-sec  Lg 
of  "sustained"  (3  or  more  cycles)  amplitudes. 

MLg  =  3.75  +  0.9  1og1QA  +  log  ]Q  (A/T)  0.5°  ^  A  £  4° 

=  3.30  +  1.66  log  ]QA  +  log1Q  (A/T)  4°  —  •A  —  30° 

Street  (1976)  and  Bollinger  (1979),  respectively,  found  Nuttli 's  formulae 
to  be  applicable  in  northeastern  and  southeastern  North  America,  provided  that 
the  maximum  distance  is  limited  to  approximately  2000  km. 

Street  et  al.  (1976),  on  the  other  hand,  assumed  C(T,  a  )  to  be  known 
and  then  specified  S(T)  such  that  the  magnitude  scales  at  different  periods 
were  set  equal  for  an  m^  =  1.5  event.  For  an  mb  =  2.5  event,  the  magnitude 
calculated  at  0.1  sec.  according  to  their  formulation  would  be  1.8,  and  the 
discrepancy  between  mb  and  -j  increased  rapidly  with  increasing  m^.  Since 
there  is  no  implicit  or  explicit  reasoning  behind  the  assumption  of  a  known 
C(T,  A  ) ,  we  are  inclined  towards  the  procedure  of  determining  C(T,  A  )  ex¬ 
perimentally  and  then  calculating  the  S(T)  so  that  a  uniform  magnitude  would 
be  obtained  at  all  periods. 


(G)  Others. 

Sn  to  Lg  conversion  appears  to  occur  near  the  margin  of  the  American  con¬ 
tinents.  For  events  from  the  West  Indies  and  Mexico  recorded  at  North 
American  stations,  Isacks  and  Stephens  (1975)  identified  the  prominent  phases 
which  arrived  after  Sn  as  possibly  a  converted  Lg  at  the  continental  margin. 
Chinn  et  al.  (1980)  observed  similar  conversions  for  events  in  the  Nazca 
Plate  recorded  at  South  American  stations.  In  neither  of  the  studies  was 
any  Lg  to  Sn  conversion  observed. 


A  number  of  investigators  have  explored  the  possibility  of  using  the 
ratio  of  Lg-ampl i tude  to  P-amplitude  as  a  discriminant  for  the  earthquake 
and  the  underground  explosion  populations.  This  possibility  was  tested  by 
Pomeroy  and  Nowak  (1979),  Pomeroy  (1980),  Nuttli  (1980  b) ,  and  Gupta  et  al . 
(1980)  for  propagation  paths  in  western  and  central  Soviet  Union,  and  by 
Pomeroy  and  Nowak  (1979)  and  Pomeroy  (1980)  for  propagation  paths  in  eastern 
and  western  United  States,  respectively.  Their  findings  indicate  a  tendency 
for  the  Lg  to  P  amplitude  ratios  to  be  greater  than  1.0  for  earthquakes  and 
less  than  1.0  for  underground  nuclear  explosions.  The  ratios,  however, 
appear  to  be  strongly  dependent  on  on  the  epicentral  distance  and  the  re¬ 
gional  attenuation  in  the  propagation  paths  and  therefore  cannot  be  used 
reliably  as  a  discriminant  between  explosions  and  earthquakes. 

Contrary  to  higher-mode  surface  waves  in  continental  structures, 
higher-mode  Love  waves  in  sediment-covered  oceanic  structures  do  not  form  a 
coherent  family  of  arrivals  at  short  periods  (Knopoff  et  al.,  1979).  This 
phenomenon  can  serve  to  explain  the  absence  of  Lg  waves  in  the  oceanic  structure. 
These  authors  also  point  out  that  since  a  large  fraction  of  the  shear  energy  at 
the  stationary  phases  of  higher-mode  Love  waves  is  concentrated  in  the  sediment¬ 
ary  layer,  absorption  by  the  low-rigidity  sediment  and  scattering  due  to  varia¬ 
tions  in  its  thickness  can  account  for  tne  rapid  attenuation  of  the  higher-mode 
Love  waves  in  oceanic  structures. 


Central  and  SE  3.65i0.04  Saint  Louis  Univ.  New  Madrid  Stauder  arid  Bollinger 

network  (Wood-  earthquakes  (1963) 

Anderson  torsion 
seisnomrter) 
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II .  Contrast  in  Lq  Wave  Propagation  in  the  Eastern  United  States  with  that  in 
the  USSR. 

The  purpose  of  this  report  is  to  portray  and  discuss  similarities  and 
differences  in  Lg  wave  propagation  in  the  eastern  United  States  with  that 
in  different  portions  of  the  USSR.  The  discussion  will  be  divided  into  five 
areas : 

1.  Lg  vs.  P  amplitudes 

2.  Lg/P  amplitude  ratios  as  a  function  of  distance 

3.  Lg  group  velocity 

4.  Lg  energy  ratios 

5.  Lg  attenuation 
1.  Lg  vs.  P  amplitudes. 

A  plot  of  Lg  vs.  P  wave  amplitudes  for  earthquakes  in  the  eastern  United 
States  is  presented  in  Figure  1.  The  dashed  line  in  this  figure  represents  a 
wave  amplitude  of  Lg  approximately  equal  to  10  times  the  wave  amplitude  of  P 
(Lg  ~  10  P).  Also  shown  in  this  figure  are: 

1.  Amplitudes  of  Lg  and  P  measured  from  records  of  the  underground  ex¬ 
plosion  SALMON  shown  by  solid  triangles. 

2.  A  solid  line  (Lg  ' ;  6.5  P)  based  on  104  measurements  of  Lg  and  P 
from  earthquakes  in  Africa. 

In  general,  this  figure  quantifies  the  general  observation  in  the  eastern 
United  States  that  Lg  is  commonly  the  largest  regional  seismic  wave  recorded 
and  often  is  the  only  signal  recorded  from  small  events. 

In  the  USSR,  the  situation  is  complicated  by  the  fact  that  most  of  the  data 
available  are  from  WWSSN  or  other  recording  stations  located  outside  the  USSR 
while  the  events  of  interest  are  within  the  USSR  usually  on  the  other  side  of 
significant  tectonic  boundaries.  With  that  in  mind,  the  data  in  Figure  2  can 
be  explained.  The  few  earthquake  events  in  the  western  portion  of  the  USSR 
shown  as  solid  squares  scatter  around  the  solid  line  represented  by  Lg  =  P/10; 
that  is,  the  amplitudes  of  Lg  are  smaller  than  amplitudes  of  P.  The  data  from 
the  Gazli  earthquakes  falls  closer  to  the  dashed  line  (Lg  =  P).  For  both  data 
sets,  the  Lg  amplitudes  are  significantly  smaller  than  those  observed  in  the 
eastern  United  States.  When  data  published  by  Soviet  authors  for  earthquakes 
in  the  eastern  USSR  is  examined,  the  situation  is  significantly  different  as 
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shown  in  Figure  3.  For  earthquakes  not  crossing  major  tectonic  boundaries 
(Events  1-12),  Lg  is  the  predominant  phase  on  the  seismogram.  Thus,  in  this 
respect,  the  eastern  portions  of  the  US  and  the  USSR  are  similar. 

For  explosions  in  the  western  central  portions  of  the  USSR,  the  data 
falls  about  or  below  the  line  Lg  =  P,  indicating  a  similarity  to  the  earth¬ 
quake  data.  This  is  shown  in  Figure  4  for  all  the  USSR  explosions  as  studied 
by  Rondout  during  this  contract.  This  indicates  little  contrast  with  the 
earthquake  population  and  suggests  that  the  propagation  path  rather  than  the 
source  properties  exert  a  predominant  effect  on  Lg  propagation  in  this  region. 
This  result  is  in  agreement  with  the  SALMON  results  in  the  eastern  United 
States. 

2.  Lg/P  ratios  as  a  function  of  distance. 

Lg/P  ratios  have  been  suggested  as  a  possible  discriminant  between 

explosion  and  earthquake  sources.  In  Figure  5,  data  on  the  logarithm  of 

the  ratio  (A/T)  Lg  are  plotted  vs.  distance  for  all  events  studied  by 
(A/T)  P 

Rondout  (with  the  exception  of  the  eastern  USSR  events).  The  earthquakes 
in  the  eastern  USSR  would  plot  as  positive  values  in  this  data.  Two  straight 
lines  represent  the  best  least-squares  fit  to  the  explosion  and  earthquake 
data.  Although  the  data  from  earthquakes  is  sparse  there  does  not  appear 
to  be  a  significant  difference  between  the  two  populations.  Because  the 
data  presented  here  comes  from  a  wide  variety  of  source-receiver  patterns, 
more  detailed  comparisons  of  the  two  populations  on  a  path-by-path  basis  are 
requi red. 

Although  similar  data  on  explosion  for  the  eastern  United  States  is 
limited  to  that  of  SALMON,  it  can  be  inferred  from  a  examination  of  Figure  1 
that  the  explosion  and  earthquake  data  would  fall  in  the  same  general  region 
on  a  plot  such  as  that  of  Figure  5. 

3.  Lg  group  velocity. 

An  initial  observation  that  the  Lg  group  velocity  for  the  SALMON  event 
was  low  led  to  a  more  extensive  study  of  the  group  velocity  values  for  Lg 
both  in  the  eastern  US  and  in  the  USSR.  Data  for  a  number  of  earthnuakes 
in  the  eastern  United  States  are  presented  in  Figure  6  (WWSSN  ana  LRSMdata) 
and  7  (NEUSSN  data).  These  show  that  the  group  velocity  from  a  number  of 
event-station  pairs  falls  below  the  generally  accepted  3.5  km/sec  value.  The 
solid  traingles  in  Figure  6  reoresent  data  from  the  SALMON  exDlosicn.  Possible 
explanations  for  these  low  group  velocities  include  the  following  two  hypotheses: 
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1.  The  propagation  paths  from  SALMON  in  general  contain  a  segment  with 

a  thick  sequence  of  sediments  in  the  Mississippi  Embayment  and  this 
thick  cover  of  low  velocity  material  may  significantly  lower  the 
average  group  velocity  of  Lg.  This  is  partially  substantiated  by 
the  low  values  of  group  velocity  observed  on  other  paths  with 
smaller  proportions  of  their  patterns  in  the  Embayment  which  exhibit 
low  values  of  group  velocity  (still  higher  than  SALMON,  however). 

2.  Because  SALMON  is  a  shallow  event,  different  relative  mode  exitation 
at  the  source  compared  to  the  deeper  earthquakes  will  occur.  The 
different  mode  excitations  could  result  in  the  waves  sampling  of 
the  near  surface  intervals  giving  rise  to  lower  group  velocities. 

To  investigate  this  possiblility  further,  data  from  earthquakes  and  ex¬ 
plosions  in  the  USSR  was  examined  and  the  results  are  presented  in  Figure  8. 

The  solid  squares  and  the  shaded  region  represent  earthquake  data  while  ex¬ 
plosion  data  are  represented  by  solid  angles.  Although  much  additional  work 
remains  to  be  done,  these  events  are  not  encouraging  for  the  use  of  Lg  group 
velocity  alone  as  a  discriminant. 

4.  Lg  Energy  ratios. 

To  investigate  the  observed  low  group  velocities  for  SALMON  further  and  to 
quantify  more  fully  differences  in  group  velocity,  an  energy  ratio  method  was 
devised.  Basically,  this  is  a  ratio  of  the  energy  arriving  in  the  oroup 
velocity  window  4.0  to  3.4  km/sec  to  the  energy  arriving  in  the  window  3.4  km/sec- 
2.8  km/sec.  The  results  from  this  analysis  indicate  that  SALMON  has  a  relatively 
low  ratio  while  earthquakes  in  the  US  have  a  higher  value.  The  energy  in  each 
group  velocity  window  was  measured  by  measuring  the  area  encompassed  by  the 
envelope  of  the  wave  train  in  a  method  similar  to  the  AR  method  as  used  by  Brune 
on  longer  period  surface  waves.  The  results  of  these  measurements  for  SALMON  and 
several  eastern  US  earthquakes  are  presented  in  Figure  9.  The  separation  between 
SALMON  and  the  earthquakes  is  clear. 

The  results  from  a  similar  analysis  on  earthquakes  and  explosions  in  the 
USSR  are  presented  in  Figure  10.  Although  the  earthquake  data  is  again  sparse, 
the  populations  seem  to  overlap  and  discriminations  is  not  achieved.  Current 
studies  involving  these  ratios  on  a  more  local  scale  may  provide  greater  under¬ 
standing  of  discrimination  capability. 
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5.  Lg  attenuation 

Lg  attenuation  has  been  the  subject  of  numerous  investigations  (see  Part  I 
of  this  report).  'Hard'  data  on  attenuation,  particularly  comparable  data  from 
different  regions,  is  still  not  readily  available.  As  part  of  the  regional  wave 
propagation  study  by  Rondout,  attenuation  measurement  in  the  eastern  US  and  in 
the  Soviet  Union  were  carried  out  and  the  results  are  presented  in  Figures  11 
and  12  respectively.  Figure  11  is  a  composite  plot  based  primarily  on  eastern 
US  earthquake  data.  Since  the  data  is  composited,  it  was  normalized  to  a  conrnon 
magnitude.  For  comparison,  Nuttli's  approximation  toy  =  0-07  deg"'*'  is  shown 
as  a  solid  line  and  the  best  least  square  fit  to  the  data  is  shown  as  the  dashed 
line. 

In  Figure  12,  Lg  attenuation  is  plotted  for  USSR  explosions.  The  ex¬ 
plosion  data  was  normalized  through  yields  assigned  by  Dahlinan  and  Israelson 
and  thus  is  subject  to  even  greater  uncertainties  than  earthquake  data.  A 

3 

straight  line  representing  an  amplitude  fall-off  of  proportion  and  1  is 
shown.  Also  shown  in  Figure  12  are  four  paths  (indicated  by  crosses)  within 
the  USSR  as  derived  by  Soviet  investigators. 

Although  a  direct  comparison  of  Figures  11  and  12  is  difficult  because 
of  different  scales,  our  tentative  conclusion  is  that,  on  the  average,  the  two 
data  sets  could  be  derived  from  the  same  population.  Work  by  Nuutli  and  Springer 
indicates  that  portions  of  the  USSR  may  have  attenuations  intermediate  between 
those  of  the  eastern  and  western  US.  Our  composites  would  average  out  these 
di fferences . 


Gneussn 
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Figure  1.  L  amplitude  vs.  P  wave  amplitude 
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Figure  2.  A  Composite  of  the  Earthquake  Data  from  Earlier  Studies  Comparing  Lg(Z)  and  P(Z) 
Amplitudes.  Note  that  the  equation  of  the  solid  line  shown  here  is  amplitude  of 
Lg  =  Amplitude  of  P/10  and  the  equation  of  the  dashed  line  is  Amplitude  of  Lg  = 
Amplitude  of  P. 
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USSR  EXPLOSIONS  AND. EARTHQUAKES. 
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Figure  12.  Lg  Attenuation  for  USSR  Explosions. 


HI.  Magnitude-Yield'  Relation  and  Others 

An  accurate  determination  of  the  magnitude-yield  relation  is  an  important 
geophysical  problem.  Aside  from  its  obvious  application  for  estimating  the  yield 
of  unknown  nuclear  tests  by  measuring  the  amplitudes  of  the  observed  seismic  waves, 
a  well-determined  magnitude-yield  relation  may  become  one  of  the  most  useful  tools 
for  calibrating  the  seismic  energy  (especially  at  short  periods)  radiated  by 
earthquakes.  The  task  of  casting  this  relation  into  a  well  defined  form,  however, 
is  not  an  easy  one.  Difficulties  can  be  traced  to  both  the  magnitude  and  the  yield 
ends  of  the  relation.  Below  we  will  describe  some  of  the  difficulties  involved. 

The  amplitudes  of  the  observed  seismic  waves  can  be  significantly  affected 
by  several  factors,  such  as  (i)  the  medium  and  the  burial  depth  of  the  source, 

(ii)  the  degree  of  seismic  coupling  between  the  source  and  the  surrounding  medium, 
and  (iii)  the  local  structures  beneath  the  source  and  the  receivers.  The  first  and 
third  factors  have  plagued  seismologists  for  years,  but  these  problems  are  cur¬ 
rently  being  solved.  To  our  knowledge,  the  second  factor  has  not  been  studied  ex¬ 
tensively,  its  effects  are  therefore  not  well  understood. 

Several  investigators  have  attempted  to  establish  the  magnitude-yield  relation 
based  on  magnitudes  that  are  determined  from  local/regional  networks  and/or  a  rel¬ 
atively  small  number  of  events.  In  view  of  the  lack  of  completeness  of  these 
studies  and  the  importance  of  this  problem,  we  have  decided  to  (i)  undertake  a 
comprehensive  compilation  of  available  published  results  that  are  relevant  to 
the  problem  of  yield-estimation,  (ii)  present  the  results  from  our  compilation  in 
a  useful  form,  and  (iii)  improve  the  determination  of  body-wave  magnitudes,  in 
a  statistical  sense,  by  increasing  the  number  of  amplitude  measurements  at  various 
epicentral  distances.  (  ISC  determines  its  body-wave  magnitudes  only  if  3  or 
more  stations  report  their  amplitudes.  It  then  applies  the  unified  magnitude  of 
Gutenberg  (1956)  to  the  amplitudes  to  determine  the  m^.  Few  stations,  however, 
have  the  habit  of  reporting  their  amplitudes  to  the  ISC.  1 


Data. 

Because  of  the  large  number  (>  400)  of  nuclear  tests  in  the  United  States 
and  the  Soviet  Union,  we  have  limited  most  of  our  data  base  to  those  underground 
nuclear  explosions  for  which  reports  on  their  estimated  yield  exist.  The  U.S. 
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data  used  is  derived  from  Springer  and  Kinnaman  (1971,  1975),  and  the  Soviet  data, 
from  Bolt  (1976)  and  Dahlman  and  Israelson  (1977).  The  magnitude  determinations 
used  are  from  Bolt  (1976)  and  the  International  Sei smological  Centre  (ISC)  Bul¬ 
letins.  There  are  some  doubts  concerning  the  source  reference  of  the  estimated 
yield  for  the  Soviet  tests,  compiled  by  Dahlman  and  Israelson,  as  well  as  the  magnitude 
of  the  Soviet  tests  as  reported  by  Bolt;  we  are  in  the  process  of  uncovering  these 
uncertainties. 

Table  IV  represents  a  compilation  of  the  U.S.  explosion  data  used  in  this 
report.  The  table  contains  the  name,  data,  origin  time,  location,  and  burial 
depth  of  the  event;  it  also  describes  the  rock-type  surrounding  the  buried  source 
(e.g.  tuff,  alluvium,  rhyolite,  etc.),  the  dimensions  (volume,  diameter,  and 
height)  of  the  collapse  cavity,  the  body-wave  magnitude  (ISC),  and  the  announced 
or  estimated  yield.  Except  for  the  magnitude,  all  the  information  was  provided 
to  Springer  and  Kinnaman  (1971,  1975)  by  the  U.S.  Atomic  Energy  Commission  (AEC). 

A  compilation  of  the  available  Soviet  data  is  outlined  in  Table  V.  This  table 
consists  of  the  date,  computed  origin  time  and  location  (Bolt,  1976),  the  body- 
wave  magnitudes  (from  ISC  and  Bolt's  compilation),  and  the  estimated  yield  for 
these  events  (Dahlman  and  Israelson,  1977). 

Based  on  the  compilations  in  Table  IV  and  V,  we  have  made  the  following 

plots: 

From  the  Soviet  data:  m^  (ISC  and  Bolt's)  vs.  estimated  yield  (Fig  13  and  14, 

respectively) 

From  the  U.S.  data:  a.  (ISC)  vs.  estimated  yield  (Fig.  15) 

b.  volume  of  collapse  vs.  estimated  yield  (Fig.  16 

c.  diameter  and  height  of  collapse  center  vs.  estimated 
yield  (Figures  17  and  18  resDecti vel y ) 

d.  volume  of  collapse  vs.  depth  of  burial  (Fig.  19) 

Information  on  the  locality  and  the  rock-type  of  the  test-site  are  also  included 
whenever  avai lable. 


Results  and  Discussion. 


A  comparison  between  the  empirically  determined  and  computed  magnitude- 
yield  relations  in  different  media  (cf.  Fig.  7-8  of  Bolt,  1976)  and  the  data 
points  in  Fig.  13  and  15  shows  that  the  U.S.  data  can  be  approximated  closely  by 
the  curve  for  granite,  whereas  the  Soviet  data  lies  roughly  between  the  curves  for 
granite  and  water.  Body-wave  magnitudes  taken  from  Bolt,  on  the  other  hand, 
show  larger  scatter  than  m^(ISC)  when  plotted  as  a  function  of  estimated  yield 
( F i g s .1 3  and  14).  There  is  some  indication  that  (i)  events  in  the  E.  Kazakhstan 
are  more  efficient  in  generating  seismic  waves  than  the  other  test  sites  of 
the  Soviet  Union,  and  (ii)  events  situated  in  tuff  and  rhyolite  generate  waves  more 
efficiently  than  those  located  in  alluvium  at  the  Nevada  Test  Site  (NTS). 

In  plotting  the  collapse  volume  vs.  the  estimated  yield  (fig.  16;,  we 
divided  the  data  into  3  groups:  the  first  two  groups  (open  and  closed  symbols) 
refer  to  events  presented  in  Figure  15,  while  the  third  group  (semi-filled  sym¬ 
bols)  consists  cf  events  that  contain  information  on  the  collapse  volume  and 
the  estimated  yield  but  not  on  the  body-wave  magnitude.  The  first  two  groups 
are  divided,  somewhat  arbitrarily,  into  normal  (closed  symbols)  and  anomalous 
(open  symbols)  events.  The  normal  events  lie  closely  together  as  a  group, 
while  the  anomalous  events  appear  to  have  unusually  small  collapse  volumes 
for  their  estimated  yields.  Fig.  17  and  18 (the  diameter  and  depth,  respectively, 
of  the  collapse  crater  vs.  estimated  yield)  were  plotted  from  the  same  data 
set.  It  is  quite  interesting  that  except  for  the  anomalous  events,  the  diameter 
of  the  collapse  crater  can  be  approximated  as  being  linearly  proportional  to  the 
logarithm  of  the  yield;  the  height  of  the  crater,  however,  appears  to  be  in¬ 
dependent  of  the  yield.  Figure  19,  which  relates  the  collapse  volume  to  the 
burial  depth,  is  composed  of  the  events  found  in  Figure  17(or  18  as  well  as 
events  without  reports  on  their  magnitude  and  yield.  This  figure  seems  to  in¬ 
dicate  three  depth-dependent  distributions:  (i)  the  volume  of  collapse  is 
independent  of  burial  depth  when  the  latter  is  less  than  about  900  ft.,  (ii) 
at  depths  between  900  and  2500  ft.,  the  logarithm  of  the  collapse  volume  is 
approximately  linearly  proportional  to  the  burial  depth,  and  (iii)  for  the 
three  events  at  deeper  than  4000  ft.,  the  volume  of  collapse  is  again  unore- 
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dictable.  A  cautionary  remark  is  deemed  necessary  at  this  point:  the  burial 
depth  of  the  test  charge  is  usually  commensurate  with  its  size;  consequently, 
the  collapse  volume  is  probably  a  complex  function  of  the  local  rock  type,  burial 
depth,  and  the  actual  yield. 

In  the  remainng  quarter  of  this  fiscal  year,  we  plan  to  conduct  three 
projects:  (i)  to  perform  regression  analysis  for  the  data  parameters  mentioned 
above,  (ii)  to  run  a  few  simple  statistical  tests  to  evaluate  the  relative  im¬ 
portance  of  the  various  parameters  on  the  amplitudes  of  seismic  waves  and  the 
dimensions  of  the  collapse  crater,  and  (iii)  to  upgrade  the  magnitude  determin¬ 
ation  by  incorporating  additional  amplitude  readings  from  stations  that  did  not 
furnish  this  information  to  the  ISC. 


Fiaure  13.  Body-wave  magnitude  (ISC)  vs.  yield  for 

_ _  —  w  a  a  r"  rti  i  rt  V* 


events  in 


line,  M.  =  0.75  log 


betweenuthese  two  q 
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Y  +  4. 345,  represents 
'n t i  ties  . 


our  preliminary. 


the  USSR.  The  dashed 
best-fitting  relation 


Y  (  kt ) 


Figure  15.  Body-wave  magnitude  (ISC)  vs.  yield  for  events  in  the  US.  Circles 
denote  tests  in  alluvium;  triangles,  tests  in  tuff;  and  rectangles,  tests  in 
rhyolite  (R),  sandstone  (Ss),  or  pillow  lava  (P.L.).  The  announced  and  esti¬ 
mated  yields  are  indicated  by  filled  and  open  symbols,  respecti vely . 


Figure  16.  Volume  of  collapsed  crater  vs.  estimated  yield  for  events  in  the  US. 
Open  and  filled  symools  refer  to  events  presented  in  Figure  3,  whereas  semi-filled 
symbols  refer  to  events  that  contain  information  on  the  collapse  volume  and  the 
estimated  yield  but  not  on  the  body-wave  magnitude  and  therefore  not  plotted  in 
Figure  3.  Filled  symbols  denote  normal  (N)  events  which  lie  closely  together  as  a 
group;open  symbols,  anomalous  (A)  events  which  aroear  to  have  unusually  small  col¬ 
lapse  volumes  for  their  estimated  yields. 


Y(kt) 


H  ( ft ): 


Figure  13.  Height  of  the  collapsed  crater  vs.  estimated  yield  for  events  in 
the  US.  Except  for  the  semi-filled  symbols,  the  letends  are  similar  to  those 
in  Figure  4. 
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Table  IV 

U.S.  Underground  Nuclear  Explosions 
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T 

9.57  E04 

500  X  30 

6 

A 

1.78  E06 

898  X  201 

300 

5.5 

A 

6.8  E04 

402  X  48 

T 

1000 

700 

6. 1 

Pillow 

9.0  E05 

2002  X  15 

1000 

6.4 

T 

330  X  20 

32 

5.6 

T 

1000  X  75 

no 

140 

5.6 

T 

17 

5.0 

T 

16 

4.9 

T 

4.7  £05 

1102  X  46 

61 

5.4 

A 

5.7  E05 

632  X  123 

30 

4.7 

T 

4.83  £05 

28 

T 

2.36  E05 

574  X  79 

20 

T 

1450  X  70 

120 

5.6 

T 

1.74  E05 

800  X  25 

25 

e 

r 

5.64  E05 

720  x  100 

25 

A 

5.91  £05 

938  X  140 

ICO 

5.3 

T 

2.7  E04 

200  X  30 

9.0 

100 

4.3 

A 

455  X  35 

6 

T 

11C0  X  65 

93 

5.5 

T 

3.54  E06 

1300  X  40 

1CG0 

1900 

6.4 

T 

6 

4.4 

T 

2.63  £05 

600  X  85 

8 

4.6 

T 

515  X  43 

6 

4.3 

T 

28 

5.0 

T 

790  X  157 

39 

T 

20 

5.1 

T 

975  X  160 

105 

no 

5.5 

T 

1010  X  175 

94 

5.5 

T 

729  X  95 

11 

T 

1100  X  100 

220 

170 

5.8 

T 

500  X  60 

10 

32 

5.1 

T 

616  X  21 

10 

T 

1000  X  73 

40 

4.9 

T 

810  X  103 

80 

100 

T 

b'a/  X  ii 

1 

Laguna 
ftatvliv  1 1 
Mlnl.it  » 


51 


U.S.  Underground  Nuclear  Explosions 


Year 

Dice  Shot  Time 

Oepth  (ft) 

Latitude  (°N) 

Longtitude 

1971 

1008 

143000 

1240 

37.11 

116.04 

cont. 

1106 

220000 

5875 

51.47 

-179.11 

1214 

210959 

1085 

37.12 

116.09 

1972 

0517 

14100 

1059 

37.12 

116.09 

0720 

171600 

1391 

37.21 

116.18 

0921 

153000 

1838 

37.03 

116.04 

0926 

143000 

970 

37.12 

116.09 

1221 

201500 

2258 

37.14 

116.08 

1973 

0308 

161000 

1866 

37.10 

116.03 

3425 

222500 

1486 

37.00 

116.03 

0426 

171500 

1850 

37.12 

116.06 

0517 

160000 

39.79 

108.37 

0605 

17C000 

1284 

37.18 

116.21 

C606 

130000 

3490 

37.24 

116.35 

0628 

191512 

1530 

37.15 

116.09 

1012 

170000 

1350 

37.20 

116.20 

1974 

0227 

170000 

37.10 

116.05 

0619 

160000 

37.20 

116.19 

0710 

160000 

37.07 

116.03 

0330 

150000 

37.15 

116.08 

0926 

150500 

37.13 

116.07 

1975 

0228 

151500 

37.11 

116.06 

C307 

15C000 

37.13 

116.08 

0405 

194500 

37.19 

116.21 

0424 

141000 

37.12 

116.09 

0514 

140000 

2510 

37.22 

116.47 

0603 

142000 

2398 

37.34 

116.52 

0603 

144000 

2090 

37.09 

116.04 

0619 

130000 

2992 

37.35 

116.32 

0 626 

123000 

4301 

37.28 

116.37 

1024 

171126 

440 

37.22 

116.19 

1023 

143COO 

4150 

37.29 

116.41 

1120 

150COO 

2630 

37.22 

116.37 

1220 

200CG0 

2349 

37.13 

116.06 

1976 

0103 

191500 

4761 

37.30 

116.33 

0204 

142300 

2100 

37.07 

116.03 

C204 

144000 

2149 

37.11 

116.04 

0212 

144500 

3999 

37.21 

116.49 

0214 

113000 

3229 

37.24 

116.42 

0309 

140000 

2851 

37.31 

116.36 

0314 

123000 

4177 

37.31 

116.47 

0317 

141500 

2984 

37.26 

116.31 

0317 

1445C0 

2559 

37.11 

116.05 

Collapse  Crater 


840  X  26 
36C0  X  55 
372  X  41 


1350  X  90 
388  X  51 
584  X  97 


Sandstone 

T 

T 

A 

T 


7 

Cathay 

6.6 

Cannikin 

24 

Chaenactis 

8 

Zinnia 

21 

4.9 

0  i  attiond 

Scu i 1 s 

no 

5.6 

Oscuro 

15 

4.1 

Gelpniniuid 

27 

4.8 

Fla* 

67 

5.3 

Mi  era 

21 

4.5 

Angus 

120 

5.6 

Slarwort 

5.1 

Hiu  bianco 

26 

5.0 

Utdo  U6cn 

570 

6.1 

Alii  el  idl'd 

60 

4.9 

Portal  ao 

9 

4.7 

Hub* /  Ace 

150 

5.6 

Latir 

20 

4.8 

Min>)  uljOe 

170 

5.7 

Escubosg 

200 

5.6 

4‘ortfc-an- 

teau 

100 

5.5 

Stan/an 

135 

5.6 

lopgul 1  ant 

120 

5.4 

Cabri 1 lo 

20 

4.9 

Dining  Car 

9 

4.5 

Edun 

380 

5.0 

T /bo 

275 

5.8 

Sti 1  ton 

160 

5.6 

Mi; Zen 

520 

5.9 

Must 

750 

6.1 

Cj.i  c-:  L«rt 

15 

4.7 

Husky  Pup 

1200 

6.2 

Kas-»eri 

500 

5.9 

Inlet 

160 

5.6 

CiiibertJ 

6n0 

6.2 

Mat- niter 

200 

5.5 

leelbOn 

150 

5.6 

Cblvdll 

900 

6.1 

fonlina 

350 

5.3 

Che  hire 

350 

5.3 

Ibt  uji'y 

900 

6.2 

Cu  !t</ 

500 

6.0 

Po«l 

200 

5.8 

Strait 

52 


Table  V 


Origin 

Year 

Oate 

Time 

Latitude  (°N) 

Longitude  (°E) 

(Bolt's) 

Announced 

Estimated 

mb  (ISC) 

1964 

0315 

75958 

49.70 

78.00 

6.2 

49 

5.6 

0516 

60058 

49.90 

78.30 

6.2 

44 

5.6 

0719 

55959 

49.90 

78.10 

6.0 

29 

5.4 

0918 

75955 

72.90 

55.20 

2 

4.2 

1025 

75959 

73.50 

53.70 

5.3 

14 

5.1 

1116 

55957 

49.70 

78.00 

6.1 

49 

5.6 

1965 

0115 

55959 

49.89 

78.97 

7.0 

125 

110 

5.8 

0303 

61457 

49.82 

78.07 

6.0 

34 

5.5 

0511 

63958 

49.79 

77.92 

6 

4.9 

0617 

24458 

49.97 

78.07 

5.8 

21 

5.2 

0917 

35958 

49.81 

78.05 

5.5 

15 

5.2 

1008 

55959 

49.89 

78.05 

5.8 

34 

5.4 

1121 

45758 

49.77 

78.06 

6.1 

47 

5.6 

1224 

45958 

49.88 

78.04 

7 

5.0 

1966 

0213 

45758 

49.82 

78.13 

6.5 

270 

6.1 

0320 

54958 

49.70 

78.00 

170 

6.0 

0421 

35758 

49.81 

78.05 

5.3 

28 

5.3 

0507 

35758 

49.74 

77.30 

4 

4.8 

0629 

65758 

49.93 

78.01 

36 

5.6 

0721 

35758 

49.70 

78.00 

5.9 

24 

5.3 

0805 

35758 

49.90 

78.00 

6.1 

29 

5.4 

0819 

35301 

50.40 

77.90 

4.6 

4 

5.1 

0907 

35158 

49.90 

78.00 

4 

4.8 

0930- 

55953 

38.80 

64.50 

5.3 

30 

1019 

35758 

49.75 

78.03 

6.3 

65 

5.6 

1027 

55758 

73.44 

54.75 

6.5 

770 

6.4 

1218 

45758 

49.93 

77.73 

6.5 

120 

5.8 

1967 

0226 

35758 

49.78 

78.12 

6.6 

210 

6.0 

0325 

55759 

49.77 

78.08 

5.9 

21 

5.3 

0420 

40758 

49.74 

78.12 

6.3 

59 

5.5 

0528 

40758 

49.81 

78.11 

6.2 

32 

5.4 

0629 

25658 

49.87 

78.10 

27 

5.3 

0715 

32657 

49.83 

78.11 

6.0 

27 

5.4 

0804 

65758 

49.82 

78.05 

5.8 

23 

5.3 

0916 

40358 

50.01 

77.82 

6.0 

22 

5.3 

0922 

50358 

50.03 

77.61 

6.0 

16 

5.2 

101 7 

50358 

49.82 

78.10 

6.1 

67 

5.6 

1021 

45953 

73.37 

54.81 

6.0 

210 

5.9 

1030 

60353 

49.34 

78.11 

6.C 

33 

5.3 

1122 

40357 

49.90 

77.30 

3 

4.8 

1208 

60357 

49.34 

78.22 

22 

5.4 

1968 

0107 

34658 

49.81 

78.02 

10 

5.1 

0424 

103557 

49.83 

78.03 

7 

5.0 

0611 

30558 

49.  S4 

78.16 

5.8 

19 

5.2 

C619 

50557 

49.96 

79.09 

6.5 

35 

5.4 

0701 

40202 

47.92 

47.95 

5.7 

46 

5.5 

0712 

120757 

49.67 

78.12 

5.9 

23 

5.3 

C320 

40558 

50.  CO 

78.00 

4 

4.8 

0905 

40557 

49.76 

78.14 

6.2 

35 

5.4 

0929 

34258 

49.77 

78.19 

6.3 

110 

5.8 

1107 

1CC205 

73.40 

54.36 

6.0 

310 

6. 1 

1109 

25358 

49.79 

78.04 

4 

4.9 

1218 

50157 

49.72 

78.06 

5.7 

14 

1969 

0307 

82658 

49.81 

78.15 

6.3 

47 

5.6 

0516 

4C257 

49.77 

78.15 

6.0 

18 

5.2 

0531 

50157 

49.98 

77.73 

C.2 

25 

5.3 

0704 

24657 

49.75 

78.19 

6.0 

99 

M 

5.2 

0723 

24658 

49.37 

78.32 

6.1 

38 

5.4 

0902 

45957 

57.41 

54.66 

5.2 

8 

11 

4.8 

0908 

45956 

57.36 

£5.11 

5.2 

8 

11 

4.c 

0926 

65956 

45.89 

42.47 

5.4 

78 

5.6 

1001 

40258 

49.81 

78.21 

5.9 

21 

5.2 

1014 

70006 

73.40 

54  61 

6.5 

340 

6.  3 

1130 

33257 

43.92 

79.00 

6.9 

160 

6.0 

1206 

70257 

43.83 

54.73 

5.7 

ICO 

5.8 

1228 

34658 

50.00 

77.82 

6.5 

72 

5.7 

1229 

40158 

49.73 

78.18 

1 

5.1 

1970 

0129 

70258 

49.80 

78.21 

5.9 

52 

5.5 

0327 

50257 

49.76 

78.31 

5.4 

1C 

5.0 

0625 

45952 

52.20 

55.69 

5.3 

5 

0623 

15753 

49.83 

78.25 

6.2 

120 

5.7 

0721 

30257 

49.95 

77.75 

6.0 

29 

5.4 

0724 

55657 

49.80 

78.17 

5.8 

21 

S.3 

0906 

40257 

49.77 

78.09 

6.0 

46 

S.H 

1014 

55957 

73.31 

55.15 

6.7 

6000 

2100 

6.6 

1104 

60257 

49.97 

77.79 

6.0 

34 

5.4 

1212 

70057 

43.85 

54.77 

6.6 

190 

6.0 

1217 

70057 

49.73 

78.13 

6.1 

35 

5.4 

1223 

70057 

43.83 

54.35 

6.6 

240 

6.0 

1971 

0322 

43258 

49.74 

78.18 

6.0 

86 

5.7 

0323 

65956 

61.29 

56.47 

5.9 

45 

51 

5.5 

0425 

33258 

49.32 

78.09 

6.4 

140 

5.3 

0606 

40257 

49.98 

77.77 

5.5 

a9 

5.5 

0619 

40358 

50.01 

77.74 

5.4 

36 

5.4 

Location 


N.Z. 

N.2. 


Uzbekistan 

N.Z. 


N.Z. 


i.  Caspian  Sea 

N  Z. 


Urals 

Urais 

I.  Caspian  Sea 
N.Z. 

.  Caspian  Sea 

Urals 

N.Z. 

E.  Caspian  Sea 
Urals 


53 


1972 


1973 


1975 


1976 


0630 

35657 

49.97 

79.05 

5.9 

0710 

165959 

64.17 

55.18 

5.1 

0919 

110007 

57.78 

41.10 

0927 

55955 

73.39 

55.10 

1009 

60257 

50.00 

77.70 

1021 

66257 

49.99 

77.65 

1022 

50000 

51.57 

54.54 

1129 

60257 

49.76 

78.13 

1215 

75259 

49.98 

77.90 

1222 

65956 

47.87 

48.22 

1230 

62058 

49.75 

78.13 

0210 

50257 

49.99 

78.89 

6.3 

0310 

45657 

49.75 

78.18 

5.8 

0328 

42157 

49.73 

78.19 

5.6 

0411 

60005 

37.37 

62.00 

4.8 

0607 

12757 

49.76 

78.17 

5.7 

0706 

10258 

49.72 

77.98 

4.8 

0709 

65958 

49.78 

35.40 

5.0 

0714 

145949 

50.00 

46.40 

3.6 

0816 

31657 

49.76 

78.15 

5.6 

0820 

25958 

49.46 

48.18 

6.3 

0826 

34657 

49.99 

77.78 

5.3 

0823 

55957 

73.34 

55.08 

0502 

85658 

49.96 

77.73 

5.3 

0921 

90001 

52.13 

51.99 

5.2 

1003 

85953 

46.85 

45.01 

6.1 

1102 

12658 

49.91 

78.84 

1124 

90008 

52.78 

51.07 

5.1 

1124 

95953 

51.84 

64.15 

5.1 

1210 

42658 

49.85 

78.18 

6.0 

1210 

42708 

50.11 

78.81 

6.7 

1228 

42713 

51.70 

79.20 

4.9 

0216 

50258 

49.83 

78.23 

5.5 

0419 

43258 

50.01 

79.72 

0710 

12658 

49.78 

78.06 

0723 

12258 

49.99 

78.85 

0315 

15958 

42.71 

67.41 

0828 

25958 

50.55 

68.39 

0912 

65954 

73.30 

55.16 

0919 

25957 

45.63 

67.85 

0927 

65958 

70.76 

53.87 

0930 

45957 

51.61 

54.58 

1025 

42658 

49.76 

78.20 

1026 

55958 

53.66 

55.38 

1027 

65957 

70.78 

54.18 

1214 

74657 

50.04 

79.01 

0130 

45658 

40  o9 

77.99 

0130 

45702 

49.83 

78.08 

0416 

55302 

•»9.99 

78.82 

0516 

30257 

49.74 

78.15 

0531 

32657 

49.95 

78.84 

0625 

35658 

49.89 

78.11 

0710 

25657 

49.79 

78.14 

0814 

145958 

58.91 

75.  SO 

0829 

95956 

73.37 

55.09 

0329 

150000 

67.23 

62.12 

0913 

30258 

49.82 

78.09 

1016 

63257 

49.97 

78.97 

1102 

45957 

70.82 

54.06 

1207 

55957 

49.91 

77.65 

1216 

62302 

49.75 

73.06 

1216 

64102 

49.82 

78.12 

1227 

54657 

49.96 

75.05 

0220 

53258 

49.82 

73.08 

0311 

54258 

49.79 

78.25 

0427 

53657 

49.99 

78.33 

0608 

32658 

49.76 

78.09 

0807 

35658 

49.81 

78.24 

0823 

35958 

73.37 

54.64 

0529 

105958 

69.59 

90.40 

1018 

85956 

70.84 

53.69 

1021 

115957 

73.35 

55.08 

1029 

44658 

49.93 

73.  )7 

1213 

45657 

49.  eo 

78.20 

1225 

51657 

50.04 

78.90 

0115 

44658 

49.87 

78.25 

0421 

50257 

49.93 

78.82 

0704 

25658 

49.91 

78.95 

0723 

23258 

43.79 

73.05 

1000 


25 

5.2 

27 

5.2 

Urals 

4 

4.5 

Urals 

770 

6.5 

N.Z. 

24 

5.3 

44 

5.5 

34 

5.2 

Urals 

34 

5.4 

3 

4.9 

210 

6.0 

N.  Caspian  Sed 

90 

5.7 

43 

5.4 

33 

5.4 

15 

5.1 

7 

34 

5.4 

Turkman 

1 

4.4 

6 

4.8 

N.  Black  Sea 

0.2 

N.  Caspian  Sea 

15 

5.0 

87 

5.7 

N.  Caspian  Sed 

35 

5.3 

690 

6.3 

N.Z. 

7 

4.9 

21 

N.  Caspian  Sed 

88 

5.6 

fU  Caspian  Sed 

350 

6.1 

11 

4.5 

Urals 

20 

5.2 

tf.  Kazakh 

70 

5.6 

620 

6.0 

3 

48 

5.5 

27 

5.4 

28 

5.2 

420 

6.1 

28 

5.3 

Uzbekistan 

14 

5.2 

U.  Kazakh 

2700 

6.8 

N.Z. 

11 

5.1 

U.  Kazakh 

210 

5.9 

N.Z. 

22 

5.2 

Urals 

19 

5.2 

7 

4.8 

Urals 

3200 

6.9 

N.Z. 

150 

5.8 

2 

4.9 

23 

5.2 

3 

4.9 

23 

5.2 

140 

5.9 

2 

4.7 

16 

5.2 

45 

5.4 

N.Z. 

870 

6.4 

•i.z. 

20 

5.0 

jrals 

15 

5.2 

43 

5.5 

1600 

6.4 

N.Z. 

2 

4.1 

3 

5.0 

6 

4.5 

51 

5.6 

77 

5.7 

30 

5.4 

60 

5.6 

35 

5.5 

14 

5.2 

550 

6.3 

N.Z. 

6 

4.8 

w.  Siberia 

1400 

6.7 

N.Z. 

700 

6.6 

N.Z. 

90 

5.8 

10 

5.1 

90 

5.7 

14 

5.2 

20 

5.3 

* 

90 

5.3 

10 

5.1 

tilla* 


i 
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